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Multiple Andreev reflections in a quantum dot coupled to superconductors: 

Eff'ects of spin-orbit coupling 
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We study the out-of-equilibrium current through a multilevel quantum dot contacted to two 
superconducting leads and in the presence of Rashba and Dresselhaus spin-orbit couplings, in the 
regime of strong dot-lead coupling. The multiple Andreev reflection (MAR) subgap peaks in the 
current voltage characteristics are found to be modified (but not suppressed) by the spin-orbit 
interaction, in a way that strongly depends on the shape of the dot confining potential. In a perfectly 
isotropic dot the MAR peaks are enhanced when the strength ur and an of Rashba and Dresselhaus 
terms are equal. When the anisotropy of the dot confining potential increases the dependence of the 
subgap structure on the spin-orbit angle 9 = arctan(a_D/Qi{) decreases. Furthermore, when an in- 
plane magnetic field is applied to a strongly anisotropic dot, the peaks of the non-linear conductance 
oscillate as a function of the magnetic field angle, and the location of the maxima and minima allows 
for a straightforward read-out of the spin-orbit angle 6. 

PACS numbers: 73.23.-b, 74.50.+r, 74.45. +c, 71.70.Ej 



I. INTRODUCTION 

Spin-orbit (SO) interaction is one of the most strik- 
ing relativistic effects that can be observed in solid state 
systems. The seminal proposal of a spin-polarized field 
effect transistor by Datta and Das^, and the range of 
possible applications in spintronics^ led to a remarkable 
progress in controlling''* and measuring^'^'^ not only the 
overall spin-orbit coupling constant^ but also the Dres- 
selhaus and Rashba ratio^i^iiiS. Recently an increasing 
interest has been devoted to studying the interplay be- 
tween spin-orbit (SO) coupling and superconductivity. 
In particular many studie a^-'^i-'^^'-'^'^i^'^'-'^^i"'^^'"'^'^ have widely 
analyzed the influence of SO interaction on the equilib- 
rium supercurrent (DC Josephson effect), pointing out 
that the multi-level nature of the dot is a crucial ingredi- 
ent. Indeed SO coupling induces, in the presence of many 
levels, inter-band spin-flip processes which can affect the 
Josephson current in a non trivial way. 

On the other hand, also non-equilibrium trans- 
port properties of nanodevices contacted to two 
superconducting electrodes have been a subject 
of intensive research since many years, from both 
theoretica l -^^i^^'^^i^^'^^i^^'^'^i^^ and experimental^ sides. 
In particular, semiconductor-superconductor hybrid 
systems are currently under spotlight, because tech- 
niques have been developed to realize highly transparent 
semi/super interfaces^, by suppressing Schottky barri- 
ers and enhancing the measured signal. In many cases 
such devices are based on In-As-related material a^^'^^ , 
which are known to exhibit Rashba and Dresselhaus 
SO interaction. Nevertheless, a systematic analysis 
of the interplay between flnite bias superconducting 
transport and spin precession induced by SO is lacking. 
The aim of the present work is to bridge this gap, 
investigating the effects of SO on electron transport 



through a quantum dot (QD) connected to two supercon- 
ducting leads, biased by a voltage V, as sketched in FiglT] 

It is well known that the subgap I-V curve of transport 
through a quantum dot is characterized by Andreev 
peaks, originating from the resonance of the dot level 
with a multiple Andreev reflection (MAR) trajectory, 
i.e. a sequence of odd Andreev reflections occurring 
at the interfaces with the leads^^. In each Andreev 
reflection the incoming and outgoing electron-hole pairs 
have opposite spins^'—. In the presence of spin-orbit, 
however, the dot levels lack of a definite spin orientation. 
The question thus arises whether the resonances with 
dot levels and Andreev processes persist, or whether 
spin-orbit coupling suppresses the MAR peaks. Another 
interesting issue is the role of the shape of the dot 
confining potential; its anisotropics indeed affect not 
only the number of effective levels involved in transport, 
but also the relative weight of Rashba and Dresselhaus 
components. The scenario is expected to be even richer 
in the presence of an in-plane magnetic field, since the 
latter introduces a preferred spin direction competing 
with the spin-orbit precession. 

Here we analyze how the subgap MAR pattern is af- 
fected by these phenomena. Our approach is based on 
the Keldysh technique, which allows us to account for 
non-equilibrium regime induced by the finite voltage. In 
order to single out the effects of spin-orbit interaction, we 
neglect Coulomb interaction in the dot, so that tunnel- 
ing between the leads and the dot can be treated exactly. 
Our results thus apply to the regime of relatively high 
dot-lead transparency (the linewidths are assumed to be 
some fraction of the superconducting gap) , so that spuri- 
ous effects due to charging are avoided. We also consider 
different shapes for the dot confinement potential, and 




V 



FIG. 1: (Color online) Scheme of the set-up under investi- 
gation: a quantum dot is coupled to two superconducting 
electrodes biased by a voltage V. 



describes the coupling with an in-plane magnetic field H 
with intensity H and angle (pn with the x axis. In Eq. (j4]) 
g is the gyromagnetic factor, fis the Bohr magneton, and 
(JxT^y the Pauli matrices. Finally 
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CTyPx) + -r- (cr^Px - CTyPy) (5) 



is the spin-orbit coupling accounting both for Rashba 
and Dresselhaus components with coupling constants 
an and ao respectively. 

We adopt as a basis the orbital eigenstates Xn (r) of the 
dot in the absence of spin-orbit coupling and magnetic 
field 



analyze the role of anisotropy. 

The structure of the paper is the following. In SecHIlwe 
describe the model. The expression for the current is pre- 
sented and briefly outlined in Sec lIIII Results are then 
presented in Sec lIVi which is divided into three subsec- 
tions, addressing the effects of spin-orbit on the subgap 
structure for an isotropic dot, the role of anisotropy and 
of an in-plane magnetic field respectively. Conclusions 
can be found in the last section V, and details of the 
calculation in the Appendices. 



II. THE MODEL 



We start by describing the isolated dot. The Hamilto- 
nian reads^ 



Ti-D = 



(1) 



drS{r) 



2m* 



Vo{r)+H„ 



l^so 



d(r) 



where d(r) = (rf|(r), d|(r))-^ is the two spin compo- 
nent electron field operator in the dot. The first term 
in the brackets of Eq. ^ describes the kinetic energy, 
p = —ihV being the momentum operator and m* the 
electron effective mass. The second term 



Vo{r) = V^^{x,y) + V^{z) 



(2) 



models the dot confining potential, for which we have 
assumed a harmonic confinement 



m 



^|(2;,y) = '-^{ijJxX^ ^ IjJyV^) 



(3) 



with frequencies LOx^y in the 2DEG x— y-plane, and a hard 
wall potential Vj^ along the growth direction z. The third 
term 



H-n 



g[iBii 

g^isH 



(4) 



{ffx cos 4>H + cTy sin 4)h) 



2m* 



V,{v)\ xn{r) = elxn{v) 



(6) 



where ej^ are the related eigenvalues, measured with re- 
spect to the equilibrium Fermi energy Ep. Typically 
transport properties within a voltage range comparable 
to |A| only involve a few dot levels. We assume that the 
width in the growth direction z is small, i.e. that the 
level spacing due to z confinement is large compared to 
|A|, so that only the lowest state along z is involved and 
that the electron dynamics is actually two-dimensional. 
Thus one has 



= EQ + h^LUxHx+UJyny) 



nx,ny = 0, 1 



(7) 



with Eq denoting the ground state level. The levels con- 
tributing to transport thus depend on the ratio of the 
confining energies hwx.y with respect to |A|. Denoting 
by N the number of these effective orbital levels, only 
2N states determine the current-voltage characteristics 
(the factor 2 arising from spin), and the electron field 
can fairly be approximated as 



d.(r) 



N 

E 



Xn(l")rf„ 



(8) 



where djj^ and dna [n = 1,...N) denote the related 
creation and annihilation operators, respectively. Insert- 
ing Eq. IH]) into Eq. H]), one obtains a 27V x 2N matrix 
representation Hd for the dot Hamiltonian 



Ti-D = 



1 "l|i ■ • ■ J "TV 



NT' U' 



d^ 



■Hr 



( di^ \ 



dii 



\dNi J 



(9) 



which can be diagonalized. 

For the leads we adopt the customary 3D s-wave BCS 

model. Labeling hy p = + (p = — ) the left (right) lead. 



one has 



k,(T=T,i 



(10) 



where cj^. denotes the electron creation operators with 
momentum k in the p-th lead and with spin projection 
a =t,ij |A| the half-gap (supposed the same in both 
electrodes), (/< the equihbrium superconducting phase 
difference between the two leads, and ^k = k^/2?7i — Ep 
the excitation energy with respect to the Fermi level Ep- 
Here Ep is the Fermi level of the two leads at equilib- 
rium; the applied bias V is included as a time-dependent 
tunneling amplitude (see below). 

Finally, a term 

N 

ntun,p = E EE(*J'^"W^k,P,A-+^-C-)(ll) 
(T=T,i k n=l 

describes the tunneling between dot and leads. We 
have assumed for simplicity that the tunneling ampli- 
tudes tp^n between the p-th lead and the level n of the 
dot are spin and k-independent. A generalization is 
straightforward. The explicit time-dependence ip.„(i) = 
^p,n(0) exp(ipwyt/2) accounts for the bias V applied be- 
tween the superconductors^^, where 



ujy 



eV 



(12) 



is the frequency related to the bias. The energies 
^p,n — 7r/9(ei?)|ip^„(0)p are the tunneling linewidths 
associated with tunneling from the n-th dot level into the 
p-the lead, with p{ep) being the DOS in the normal state. 

In conclusion the total Hamiltonian of the system reads 

'H = Hd + 'Hlead,+ + 'Hlead,- + 'Htun,+ + 'Htun,- (13) 

where Hd-, 'Hiead,±, and Titun,p are given by Eqs.Q, ^Q 
and (jlip respectively. 

III. CURRENT 

The current flowing from p-th lead into the dot can be 
written as 

^p(^) = ^ E E E (.tpAt){cl.d„.) - h.c.) (14) 

a=T,i k n=l 

Combining the non-equilibrium Green function technique 
with the Dyson equation, one can rewrite Eq. (I14p as 

Ipit) = -2e^ f dt'TT4N {{a,)AN pp(t, i')G(i', t)]+"](15) 



Here bold notations denote matrices in the Keldysh space 
of the quantum dot, with the superscript -l-(— ) labeling 
the upper (lower) Keldysh time contour branch. In par- 
ticular Up and G respectively describe the self-energy due 
to p-th lead and the dot Green function, evaluated in the 
presence of the leads 



G-(Go 



Z_) 



(16) 



where Go describes the isolated dot. Details of the 
derivation of Eq. (|15p can be found in ApplXl Here 
we emphasize that, differently from the customary 
treatment of superconducting transport through a 
quantum dot, in our case one has to adopt a Nambu 
space for the dot with dimension 4N, instead of 2N, 
due to the spin-orbit term which effectively induces 
spin-flip tunneling processes. For these reasons in 
Eq. p5|) the symbol Tr47v denotes the trace on the 
4Ar-dimensional (extended) Nambu space of the dot, 
and (CTz)4Ar = (T^ (g) l2Ar, where az is the usual Pauli 
matrix and I2N the 2A^ x 2N identity matrix. 

The presence of a finite bias yields the current to be 
time-dependent (AC Josephson effect); it is thus suit- 
able to derive a harmonic time series representation of it. 
By applying a discrete Fourier Transform (see App|B|) to 
llp{t,t') and G{t',t), one obtains 



00 
Ip{t)^~-n ^ 

m— — oo 



^-i-imLUvt 






dLoTliNicTz) 



AN 



.ni,n2 — — 00 



where 



IIp(ni, 712; uj) G(n2, rii -t- m; uj) 






(17) 



(18) 



is the voltage-dependent fundamental domain, and m la- 
bels the harmonic. The Fourier coefficients lLp{ni, 712; lo) 
and G(ni,n2;w) of the self energy and the dot Green 
function are computed in App[Bl We emphasize that, 
since Eq. (|17[) is formally exact, charge conservation im- 
plies that 



/+(t)--/_(t)-/(t) 



(19) 



In the explicit calculation the series in ni and n2 must 
be truncated to a cut-off, which is chosen in order to 
ensure convergence of I± and the fulfillment of Eq. (|19p . 

The DC component of the current P^ corresponds to 
the m — Q harmonic and thus reads 



-p-5R / dwTr4Ar(cr^)4Ar 
TT Jf 



E 

_Tli ,712 — — 00 



Tp{ni,n2;Lo)Q{n2,ni;uj) 



(20) 
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FIG. 2: (Color online) Sketch of the dot confining potential 
Q for different values of the frequencies in the x and y direc- 
tions, a) huix,y 3> |A| : only one level takes part to transport; 
b) hux = huiy ~ |A|: multilevel isotropic dot; c) huix — |A| 
and hujy 3> |A|: multilevel strongly anisotropic dot. 



where p = ±1. While the case of equilibrium Josephson 
current {V = and 7^ 0) was thoroughly analyzed in 
Reflla. here we focus on the out of equilibrium case {V ^ 
0), and henceforth set = 0. The current (PO)) is then 
evaluated with numerical integration in the frequency. 



z-axis, because the Dresselhaus term in Eq. ([5]) does not 
commute with the total angular momentum J^ — Lz+Sz- 

We start by some general remarks. In the first 
instance, simple perturbation theory arguments lead to 
conclude that spin-orbit effects can only be observed in 
a multilevel quantum dot: for one single level quantum 
dot, spin-orbit interaction does not play any role, be- 
cause the momentum operator appearing in Eq. ^ can 
only couple dot levels with different quantum numbers 
nx and Uy (see Eq. ([7]))^. For these reasons, in the 
following we shall consider only multilevel quantum 
dots. The number 27V of states contributing to transport 
depends on the parameters of V\\{x,y), and is chosen in 
such a way that the subgap I-V curve remains unaffected 
by inclusion of higher energy levels. Notice that for 
intermediate level spacing (i.e. for huj^.y > |A|) the dot 
excited levels are energetically too high to contribute to 
transport directly, but close enough to the unperturbed 
ground state to modify its resonance conditions via spin- 
orbit interaction, hence affecting transport indirectly. In 
this sense the effect of spin orbit can be considered as 
due to a single level. 

Secondly, the shape of the confining potential has a 
crucial role in determining spin-orbit effects: it is easy to 
verify that entries of the dot Hamiltonian ([1]) originating 
from spin-orbit couplings ([5]) scale as olr^d /^x.yi where 



At 



rWu^^y 



(21) 



are the lengths associated with the frequencies of the 
confining potential Vn. Indeed the value of uox/^y 
determines the number of states in the x and y direction 
that contribute to transport, and also changes the 
relative weight of spin-orbit terms proportional to px 
with respect to those involving py. 



IV. SPIN-ORBIT AND MULTIPLE ANDREEV 
REFLECTIONS 

Here we present the results concerning the DC 
component (PO)) of the current, analyzing the effects of 
spin-orbit interaction on the multiple Andreev reflections 
pattern and subgap structure. We consider different 
types of shapes for the confining potential V\\, from 
isotropic {lOx = ^y) to strongly anisotropic (10 y 3> Wx), 
as shown in FigO Experimentally, a tuning of the 
frequencies ujx,y (see Eq. ([3])) can be achieved through 
side gates, which do not substantially alter Rashba and 
Dresselhaus spin-orbit coupling. For the sake of clarity, 
it is worth pointing out that throughout the paper we use 
the term isotropy as referring to the confining potential 
V|| only, and not to the full dot Hamiltonian. It is indeed 
easily checked that, even for a isotropic potential V\\ 
and in absence of in-plane magnetic field, the Hamilto- 
nian of the dot is not invariant under rotation around 



Finally, since we account for both Rashba and Dres- 
selhaus terms, it is worth introducing a total Dressel- 
haus/Rashba coupling constanti^ 



'■D 



'R 



(22) 



and a relative Dresselhaus/Rashba spin-orbit angle 6, de- 
fined through the relations 

cost^ = — smfc' = . 

a a 



(23) 



With the help of Eqs.(122]) and (|23l), the spin-orbit cou- 
pling ([5]) can be rewritten as 



a ^ 

Vso == -rP\\ 



(24) 



where pu — (jJx,Py) is the x — y plane momentum and 
a = (flx, ciy) is a vector with components 



ax 


= (7x sin 9 — ay cos 6 


(25) 


ay 


= (Jx cos 9 — <7y sin 9 


(26) 




FIG. 3: (Color online) Pictorial scheme of multiple Andreev 
reflections in a multilevel quantum dot. The left and right 
parts represent the energy spectrum of the superconducting 
leads, characterized by a gap |A|, whereas the horizontal lines 
in the central part of the figure describe the dot levels. The 
subgap MAR processes related to eV^/2|Aj = 1/n with n — 
3, 5, 7 are depicted. 



Eq. ([24)1 describes the well known effect that, in the 
presence of spin-orbit couphng, the spin orientation 
depends on the momentum p|| direction, in a manner 
which is also related Rashba/Dresselhaus angle 6. For 
the case of a two-dimensional electron gas, this aspect 
was recently emphasized e.g. in RefiOJ. In a quantum 
dot, however, states with well defined momentum are 
not eigenstates, due to the presence of the confining 
potential Vh . Hence, in the presence of spin-orbit 
coupling the dot levels do not exhibit a definite spin 
orientation in general, and the question arises whether 
MAR resonances are suppressed by spin-orbit. 

In the rest of the paper we present our results address- 
ing the following questions: i) does spin-orbit interaction 
in a quantum dot suppress MAR peaks? (see Figl3]) ii) 
does the geometry of the dot matter in observing spin- 
orbit effect on electron transport; iii) what is the effect 
of the interplay between the spin-orbit angle 9 and the 
magnetic field on the I-V curves? 



A. Spin-orbit effects for an isotropic dot 



We start our analysis from the case of an isotropic 
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FIG. 4; (Color online) Current (units of 2e|A|//i) as a func- 
tion of the source-drain bias for a three-level quantum dot 
with isotropic confining potential (fiwi, = hwy — 20|A|) and 
linewidth T±^i = |A|/4 (i — 1,2,3). The solid curve refers 
to the case of spin-orbit coupling a = 10~^^eVm (7 = 0.174) 
and spin-orbit angle 6 = 7r/4. The dashed curve describes 
the case of vanishing spin-orbit coupling a = 0. The spin- 
orbit coupling does not suppress the MAR subgap pattern; it 
modifies the location and the number of the peaks. 



a dot with £^0 = and confining potential frequencies 
huJx = fi^y = 20|A|, so that only the ground state lies 
in the gap energy window. The solid curve is obtained 
in the presence of spin-orbit coupling a — 10~^"'^eVm 
and Rashba/Dresselhaus angle 9 — 7r/4. For the sake 
of comparison we have also plotted the case of vanishing 
spin-orbit (dashed curve) . Our result indicates that spin- 
orbit coupling modifies, rather than suppress, the MAR 
pattern. This is due to the fact that, although spin-orbit 
coupling does not allow electrons in the dot to have a 
definite spin orientation, it does not break time-reversal 
symmetry (TRS). As a consequence of TRS, the dot lev- 
els are always doubly degenerate and the related pair 
of eigenstates mutually connected by the time-reversal 
transformation 



T: 



iayK 



(27) 



confining potential (cu^ 



CO. 



yi^ 



without magnetic field. 



The solid curve of Fig[5] shows the I-V characteristics for 



where K is complex conjugation. An electron emerging 
from an Andreev reflection at the interface can thus tun- 
nel into the doubly degenerate dot level, independent of 
its original spin-direction, even in the presence of a finite 
spin orbit coupling a. With respect to the case a = 0, the 
electron simply redistributes differently between the two 
TRS dot states. Although spin-orbit effect does not de- 
stroy the MAR subgap pattern, it does modify it, because 
the location of MAR peaks and their sharpness depend 
on the dot resonant levels and the dot-lead linewidths, 
which are both affected by spin-orbit coupling. The sim- 
ple case of an isotropic quantum dot with 3 levels (one s 
and two orbitally degenerate p levels) can elucidate this 



effect, since the spectrum can be evaluated analytically, 
obtaining 



£i — Ef) 
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(29) 
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h V hujx 
is the dimensionless spin-orbit coupling. 



(31) 



In conclusion, the MAR pattern is quantitatively 
modified by spin orbit coupling a, although it quali- 
tatively resembles the one of a dot with appropriately 
renormalized levels and linewidths. 

Let us now discuss the effects of the spin-orbit angle 6. 
In FigO the current-voltage characteristics is shown for 
different values of 9. One can notice that for an equal 
weight of Rashba and Dresselhaus terms {6 = 7r/4) the 
MAR peaks are enhanced with respect to the two cases of 
purely Rashba {6 = 0) and purely Dresselhaus {6 — 7r/2) 
interactions. Indeed for the particular value 9 — 7r/4 
electrons in the dot have a well defined spin orientation, 
since the Hamiltonian ([l} commutes with the spin oper- 
ator (ax — (Ty)/v2- In contrast, for arbitrary values of 
9 spin is not a good quantum number. The curves at 
9 = and 9 = tt/2 turn out to coincide, since the Hamil- 
tonian corresponding to these two cases can be mapped 
into each other by the unitary transformation on the spin 
variables 



V 



(32) 



B. 



The effect of the anisotropy in the confining 
potential in the absence of magnetic field 



The equivalence discussed above between transport 
properties of a dot with purely Rashba and purely Dres- 
selhaus terms persists also in the presence of anisotropy 
in the confining potential Vj| of the dot. However, when 
the anisotropy of the dot increases, the dependence of the 
I-V curves on the angle 9 becomes weaker (see FiglHl). In 
particular, for a dot with a strongly anisotropic confin- 
ing potential (e.g. for fkUy ^ huJx, |^|)j and in absence 
of magnetic field, electron transport becomes insensitive 
to the spin-orbit angle 9, and only the total intensity a 
of spin-orbit coupling matters. In order to understand 
this effect one can observe that, due to the anisotropy of 
the confining potential, only the lowest quantum num- 
ber related to the y effectively matters. Thus the opera- 
tor py appearing in the spin-orbit interaction in Eq. ([S]) 
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FIG. 5: (Color online) Current-voltage characteristics for a 3- 
level quantum dot with isotropic confining potential (fitJj; = 
hujy = 4|A|) and linewidths r±,i = |A| {i = 1,2,3), in 
the presence of a spin-orbit coupling constant 7 = 0.5 (see 
Eq. (|3ip ) and for different values of the Rashba/Dresselhaus 
angle 9. The solid curve refers to 6 = ■7r/4, i.e. equal weight 
of Rashba and Dresselhaus terms, whereas the dashed (dot- 
ted) curve refers to 9 = {9 = 7i"/2), corresponding to purely 
Rashba (purely Dresselhaus) coupling. The mixed coupling 
enhances the MAR peaks, whereas the two other cases turn 
out to coincide, due to symmetry reasons (see text). 



does not alter the dot levels involved in transport, and 
can effectively be dropped from the Hamiltonian. Impor- 
tantly, this implies that for a strongly anisotropic quan- 
tum dot the spin-orbit interaction does not induce any 
spin-precession, but rather determines one preferred ori- 
entation of the electron spin in the x — y-plane. Although 
such spin direction formally depends on the spin-orbit an- 
gle 9, the latter can indeed be gauged away completely 
from the Hamiltonian by performing a unitary transfor- 
mation only on the spin degrees of freedom 



V« 



exp 



ii--9)aj2 



yielding a spin-orbit coupling 

a 



Vf 



so 



-PxCTx 



(33) 



(34) 



which involves only the x-component of spin. In absence 
of magnetic field all directions are thus equivalent for 
electron transport, and this explains the lack of depen- 
dence on 9 of the I-V curves. 



C. Effects of the in-plane magnetic field 

In this section we discuss the effect of the in-plane mag- 
netic field, analyzing both the dependence on its inten- 
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FIG. 6: (Color online) Current (units of 2e\A\/h) as a 
function of the source-drain bias, for a 3-level dot in an 
anisotropic confining potential {ujy = I6CJ3;), with dot-lead 
hnewidths r±,i — |A| (i — 1,2,3) and spin-orbit coupling 
constant 7 — 0.5 (see Eq. (JSTJ), for different values of 
the Rashba/Dresselhaus angle 6. The solid curve refers to 
9 = it/4:, i.e. equal weight of Rashba and Dresselhaus terms, 
whereas the dashed (dotted) curve refers to 6 — {9 = 7r/2), 
corresponding to purely Rashba (purely Dresselhaus) cou- 
pling. With respect to the isotropic case, the dependence on 
9 is much weaker in the presence of an anisotropic confining 
potential. 
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FIG. 7: (Color online) Efi^ects of the magnetic field on the 
MAR subgap structure. The case of an isotropic 3-level dot 
with _Bo = 0, level spacing huj^.y ~ 4|A|, spin-orbit parame- 
ters 7 = 0.5 and 9 — it/4:, and dot-lead linewidths T±^i = |A| 
(j = 1, 2, 3). The solid curve describes the case without mag- 
netic field {B = 0), whereas the dashed and dotted curves 
refer to an applied magnetic field along the x-axis with in- 
tensity B — 0.5 and B — 1 respectively. The presence of 
a magnetic field suppresses the current at low bias and the 
MAR peaks. 



sity H and on the angle (J)h ■ FiglT] shows the subgap 
structure for different values of the dimensionless mag- 
netic field 



B 



gfiBH 



(35) 



in an isotropic 3-level quantum dot. As one can see, 
MAR peaks are suppressed by increasing B, as expected 
since magnetic field breaks time reversal symmetry and 
suppresses the probability that an electron emerging from 
an Andreev reflection matches resonance conditions with 
the appropriate spin direction. 

A more interesting scenario emerges when analyzing 
the dependence of the subgap structure on the magnetic 
field angle (pn- To this purpose, it is worth discussing 
how the spectrum of the dot varies with (/>// while keep- 
ing the intensity H constant. The shape of the confining 
potential turns out to play a crucial role on it. We start 
with the case of an isotropic quantum dot with three or- 
bital levels (one non-degenerate s level and one doubly 
degenerate p level). Figl5^ shows (some of) the levels of 
a dot with level spacing fkOx — fnoy = |A|/2, in the pres- 
ence of a spin-orbit couphng 7 = 0.5 (see Eq. (PT|) ). and 
under a magnetic field of intensity B = 0.5 (see Eq. (|35l) ). 
As one can see, the eigenvalues oscillate as a function of 



(j)H and exhibit maxima and minima located at (J)h = 7r/4 
and (l)H — 37r/4. Notice that in this case the spin-orbit 
angle 6 only affects the amplitude of the oscillations and 
not the location of the minima and maxima as a function 
of (pH- The origin of this effect boils down to the sym- 
metry of the Hamitonian. Indeed for an isotropic dot the 
Hamiltonian characterized by a magnetic field with an 
angle (pH = t^/'^ ~ 54>h can be mapped into the one with 
an angle (pn = 7r/4 -|- 5(J)h through the unitary transfor- 
mation 
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(36) 



and it therefore exhibits the same spectrum as the latter. 
For an isotropic dot, as far as the dependence of I-V 
curves on the in-plane magnetic field is concerned, the 
role of is not qualitatively different from the one of 
spin-orbit strength a, and it is thus difficult to extract 
information about spin-orbit angle by the analysis of the 
MAR pattern as a function of (pn ■ 

The situation is quite different for a strongly 
anisotropic confining potential V\\ . The dependence of 
(some) dot levels is shown in FiglSj) for a 3-level dot with 
spin orbit coupling 7 = 0.5 and magnetic field intensity 
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FIG. 8: (Color online) Behavior of (some) dot energy levels 
as a function of the angle (pH of the in-plane magnetic field 
with dimensionless intensity B — 0.5 (see Eq. (|35|l '). and for 
different values of the spin-orbit Rashba/Dresselhaus angle 8 
(solid curves refer to 6 = 7r/4 and dashed curves to 6 = 7r/20). 
a) case of an isotropic 3-level dot with hcux = hcuy — | A|/2; the 
maxima and minima of the dot level oscillations are located 
at (f>}j = 7r(2n + l)/4 (n = 0, 1,2,3) and are independent of 
the spin-orbit angle 6, which determines only the oscillation 
amplitudes, b) case of a strongly anisotropic 3-level dot with 
hujx — |A|/2 and ftLi;y/|A| — > oo. The maxima and minima 
of the dot level oscillations are located at the spin-orbit angle 
(f>H = 0±mi/2 {n = 0, 1, 2, 3). Vertical dotted lines are guide 
for the eye. 



B = 0.5. As one can see, the maxima and minima do de- 
pend on the spin-orbit angle, and are precisely located at 
^^ = 9±mr/2 (n = 0, 1, 2, 3). To understand this effect, 
we recall that when the anisotropy of the dot confining 
potential is strong, the electron exhibits a well defined 
spin orientation which depends on 9; the introduction of 
a magnetic field simply induces a spin precession around 
its direction'^^. More formally, the transformation 
yields an effectively rotated magnetic field 



H 



H{cos{ 



>H 



n/2)ax 



sm{(j}H -0 + 7r/2)cry) 
(37) 
whose direction depends on the spin-orbit angle 6. The 
structure of Eq. (|37p shows that, in the presence of an 
in-plane magnetic field, there are resonances (antireso- 
nances) whenever the angle 4)h of the physical magnetic 
field equals 0— 7r/2 (0), corresponding to an effective mag- 
netic field pointing in the parallel (orthogonal) direction 
as the spin-orbit term ([M|) . A simple example can illus- 
trate this effect: let us consider the case where only two 
levels in the x direction are present. In this case the 4 
non-degenerate eigenvalues of the dot can be computed 
analytically and read 



(38) 
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FIG. 9: (Color online) Non-linear conductance G — dl /dV 
(in units of 2e} jK) as a function of the applied bias V (units 
of |A|/e) and of the angle (\)h of the in-plane magnetic field 
with intensity B = 0.7 for a strongly anisotropic quantum dot 
{Jvjj^ = |A| and fujjyl\^\ -^ co) with two levels, and dot-lead 
linewidths r±,i,2 = |A|/2. a) without spin-orbit coupling; 
b) in the presence of a spin orbit coupling with parameters 
7 = 0.3 and angle d = 0.37r. The spin-orbit coupling yields os- 
cillations of the conductance in </>_fr with minima and maxima 
given by the spin orbit angle Q and Q — 7r/2. 



where 7 is the dimensionless spin-orbit coupling (|3ip . 
and h = g^.BH/fiuj. The resonances between 9 and (pn 
are described by the last oscillatory term in Eq. (|38p . 

A novel effect thus arises for a strongly anisotropic dot 
under an in-plane magnetic field: by analyzing the non- 
linear conductance G = dl/dV as a function of the bias V 
and the magnetic field angle (jj^i one can see that in the 
presence of spin-orbit G exhibits maxima and minima 
located at 4'h — 9 and (pH = 9 — tt/2, i.e. when the 
magnetic field direction matches the Rashba/Dresselhaus 
spin-orbit angle. The analysis of the oscillations of the 
MAR peaks as a function of (J>h allows to gain the spin- 
orbit angle 9, as illustrated by the contour plots in FiglH] 
for quantum dot with two levels. The upper panel refers 
to the case without spin-orbit coupling, where no oscilla- 
tions as a function of (J)h are present, since the magnetic 
field simply determines the preferred spin-direction. The 
lower panel describes the case of SO coupling 7 = 0.3 
and 9 = 0.3tt. In this case the magnetic field direction 
interplays with the SO angle, so that maxima and min- 
ima appear at (Jjh = ^ and (J)h = 9 — 7r/2, as indicated 



on the right of the figure. 



APPENDIX A: EVALUATION OF THE 
CURRENT 



CONCLUSIONS 



To conclude, we have investigated the effect of Rashba 
and Dresselhaus spin-orbit couphngs on the out-of- 
equihbriuni transport properties of a quantum dot cou- 
pled to two superconductors. We have analyzed how the 
I — V is affected by the total coupling constant a (see 
Eq. (P^ ) and by the Rashba/Dresselhaus angle 9 (see 
Eq. (|23|) '). discussing also the role of the anisotropy of 
the confining potential and the in-plane magnetic field 
for multi-level dots. We have found that, although SO 
effect prevents electrons tunneling in the dot to have a 
definite spin-orientation, the MAR subgap structure is 
not suppressed by SO interaction, due to the conservation 
of time-reversal symmetry. The MAR pattern is never- 
theless quantitatively modified by a, which changes the 
resonance conditions and the linewidths, affecting the lo- 
cation and the number of the Andreev peaks (see FigU]). 
The role of 9 is strongly dependent on the shape of the 
dot confining potential and the magnetic field. In absence 
of magnetic field, an isotropic dot with equal weights for 
Rashba and Dresselhaus terms exhibits higher MAR cur- 
rent peaks than the cases with purely Rashba and purely 
Dresselhaus terms (see Fig 15]), whose I-V curves coincide 
due to Hamiltonian symmetry. By increasing the dot 
anisotropy the dependence on 9 vanishes, and the MAR 
peaks are only affected by the total coupling constant a 
(see Figl6]). When an in-plane magnetic field is applied, 
the scenario is even richer, especially as a function of the 
angle (pH between the magnetic field direction and the 
X axis. For an isotropic potential the dot energy levels 
oscillate as a function of 0//, with maxima and minima 
located at (j)*^^ = 7r(2n-|- l)/4 (n = 0, 1, 2, 3), independent 
of the spin-orbit angle 9; the latter only determines the 
oscillation amplitudes and plays a similar role as the cou- 
pling a. In contrast, in a dot with strongly anisotropic 
confining potential the maxima and minima are located 
at magnetic field angles 0^ ~ 9 ± mr/2 [n — 0, 1, 2, 3). 
This enables a direct read-out of the spin-orbit angle 
through the inspection of the non-linear conductance as 
a function of the bias and the angle of the magnetic field 
(l)H (see FigE]). 



In this appendix we show how to obtain Eq. (JTSJ) from 
Eq. pi)) . As a first step, we observe that, due to the spin- 
orbit term, states of different orbital levels and opposite 
spins are coupled, and the Hamiltonian matrix ([9]) can 
be diagonalized through a unitary transformation U on 
the dot level operators 



d„,^{U,)n,D^ j^l,...2N (Al) 



obtaining 



2JV 



Hd^Y. ^Md, 



(A2) 



where Si {i — 1,2N) are the eigenvalues. In Eq. (|A1|) 
U| (Uj) is the N x 2A^ upper (lower) submatrix of the 
unitary transformation U. Introducing a block-diagonal 
2 X AN tunneling matrix 



[tp^i, ■ . . ^tppfjU-f 







(A3) 



with time-dependence given by 

/ ipuJvt/2 n \ . 

Tp{t)^^^ Q ^-.pL/2J-rp(0) (A4) 

one can rewrite Eq. (fTTj) in a Nambu notation as follows 



Ti-tun.p — 7 , ( '^k.p.f'^-k.p.l ) •'^p(^) 



£>. 



2N 



D] 



+ h.c. 



Similarly, the current (I14p can be rewritten as 



Ip{t) 



(A5) 



^5R^Tr4Ar{(a.)4iv [Tt(t)Gkp,o(i,<)]^ } 
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where 



and 



Tt ^I^P 

p ^ -tpt 



C _ / ^kp.O ^kp.O 

'-'kp.O '-'kp^O 



(A6) 



(A7) 
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denotes the lead-dot Green functions with entries 



^G;:;,:^(ii,i2) = 



CkpT 



{t,){D^['^\-,Dif,D[ 



V2) ri^V2] 



(A8) 



ih) 



with 771^2 = ± Kcldysh labels. Exploiting Keldysh-Dyson 
equation 

1 r°° 

Gkp,o(i, ^) = )^ / gkp,kp(i, i')Tkp(t') G(t', i) dt' (A9) 

one can express Gkp,o in terms of the Green function 
gkp.kp of the p-lead and the Green function G of the dot. 
The latter are defined similarly to Eq. (|A7[) with entries 



^9'^Xiti,h) 



and 



( c'"^^ 

' '^kpT 
\ '^ -kpi 



(AlO) 
(All) 



D 



2N 



D] 



ih) {dI ... dIj^ D, ... D^j,) ih) 



\DlJ 



respectively. Inserting Eq. (|A9|) into Eq. (|A5|) and recall- 
ing that the self-energy due to the p-th lead reads 

^p{ti,t2) = ^J2 Tkp(^l) gkp,kp(tl, i2)Tkp(i2) (A12) 



APPENDIX B: DISCRETE FOURIER 
TRANSFORM 



In this appendix we sketch the calculation that leads to 
the evaluation of llp(ni,n2;uj) and G(ni,n2;uj) appear- 
ing in the current formula (J17p . One first introduces^ 
a discrete Fourier transform f{n,m;u}) of an arbitrary 
two time arguments function f{ti,t2), defined through 
the relation 



f{tl,t2)-- 

E 



2^ 



(Bl) 

-i{u+nuv)ti ^i(uj+mujv)t2 f f^ m,' io) 



n.m — — C30 

where F is the fundamental domain introduced in 



Eq. p^ . The inversion formula reads 



2tt5{uji - UJ2) /(«!, 712; wi) 



(B2) 



dt 



dt' fit t') g''('^i+"i'^V')* -i{uJ2+n2U}v)t' 



where loi,uj2 G F. For a function that depends only on 
time differences /(ii, ^2) = /(^i — ^2), the discrete Fourier 
transform reads 

f{ni,n2]uj) = 5n^^n2 K'^ + niujv) (B3) 

where /(w) is the usual Fourier Transform. 

The function lLp{ni^n2]Lo) is thus obtained by insert- 
ing Eq. (|XT2l) as / into the right-hand side of Eq. (|B2|) . 
and making use of Eq. (jA4p and of the definition 



,(w 



/+00 
rfie^"*Egkp,kp(t) (B4) 

-00 1. 



one obtains Eq. p5 



obtaining 



5n2,nt[gp{uJ + riiUJv - ^)]ll 5n2,ni-p[gp{i^ + niUv - ^)]l2 

lLp{m,n2-.Lo) = ^Tt(O) I |Tp(0) (B5) 

<5„2,ni+p[gp(w -I- niwy 4- ^)]21 (^n^,,!! [gp(cj + niuv + ^)]22 



1 

fi2 



where [. . .].y (i, j = 1,2) denote the entries in the 2-dimensional Nambu space of the leads. A standard calculation 
yields'^^ 



gp{uj) ^ TTpiep) 






1 _^ 

huj 



(B6) 



where fp denotes the Fermi function of the p-the lead, and 



hi{uj) = 



-hLue{\Ap\ - \io\) 

x/|ApP-a;2 



(B7) 
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h2iuj) 



IHulQi 



V^^-|ApP 



(B8) 



with O being the Heaviside function. In Eq. (jB6[) the 
symbols (. . .)n and , (. . .)k denote matrices acting on 
the lead Nambu and Keldysh space respectively. 



means of the second Dyson equation ((T6| , where Go de- 
scribes the green function of the isolated dot. Explicitly 

GQ^{ni,n2;uj) ^ 5„^^n2'^o^i^ + niujv), where 

G^^(w) = wl4Ar-cr^(g)diag(ei,...,£2Jv) (B9) 



Finally the dot Green function G can be evaluated by with 14^? the AN x AN identity matrix. 
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